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An  estimate  of  the  Kelvin  impulse  of 
a  transient  cavity 

By  J.  P.  BEST'  and  J.  R.  BLAKE* 

1  Materials  Research  Laboratory  (MRL)  -  DSTO.  PO  Box  50.  Ascot  Vale. 

Victoria.  3032,  Australia 

1  School  of  Mathematics  and  Statistics.  University  of  Birmingham.  Edgbaston. 

Birmingham.  BI5  2TT,  UK 

(Received  9  April  1992  and  in  revised  form  18  August  1993) 

The  Lagally  theorem  is  used  to  obtain  an  expression  for  the  Bjerknes  force  acting  on 
a  bubble  in  terms  of  the  singularities  of  the  fluid  velocity  potential,  defined  within  the 
bubble  by  analytic  continuation.  This  expression  is  applied  to  transient  cavity  collapse 
in  the  neighbourhood  of  boundaries,  allowing  analytical  estimates  to  be  made  of  the 
Kelvin  impulse  of  the  cavity.  The  known  result  for  collapse  near  a  horizontal  rigid 
boundary  is  recovered,  and  the  Kelvin  impulse  of  a  cavity  collapsing  in  the 
neighbourhood  of  a  submerged  and  partially  submerged  sphere  is  estimated.  A 
numerical  method  is  developed  to  deal  with  more  general  body  shapes  and  in 
particular,  bodies  of  revolution.  Noting  that  the  direction  of  the  impulse  at  the  end  of 
the  collapse  phase  generally  indicates  the  direction  of  the  liquid  jet  that  may  form,  the 
behaviour  of  transient  cavities  in  these  geometries  is  predicted.  In  these  examples  the 
concept  of  a  zone  of  attraction  is  introduced.  This  is  a  region  around  the  body,  within 
which  the  Kelvin  impulse  at  the  time  of  collapse,  and  consequent  jet  formation,  is 
expected  to  be  directed  towards  the  body.  Outside  this  zone  the  converse  is  true. 


1.  Introduction 

Upon  the  collapse  of  a  transient  cavity  in  the  neighbourhood  of  boundaries,  a  liquid 
jet  is  often  observed  to  form  and  completely  penetrate  the  bubble  (Benjamin  &  Ellis 
1966;  Lauterbom  &  Bolle  1975;  Vogel,  Lauterbom  &  Timm  1989).  Computations  of 
the  bubble  shape  are  in  agreement  with  these  observations  (Plesset  &  Chapman  1971 ; 
Guerri,  Lucca  &  Prosperetti  1981 ;  Blake,  Taib  &  Doherty  1986).  Significant  attention 
has  been  paid  to  these  jets  as  their  impact  against  solid  boundaries  is  a  cause  of 
cavitation  damage.  It  is  also  evident  that  the  pressure  gradient  induced  by  a 
gravitational  field  is  a  cause  of  jet  formation  such  as,  for  example,  occurs  in  the  case 
of  the  bubble  formed  by  an  underwater  explosion. 

The  concept  of  the  Kelvin  impulse  was  invoked  by  Benjamin  &  Ellis  (1966)  in  order 
to  explain  the  deformation  from  a  spherical  shape  of  a  translating  cavity,  the  formation 
of  jets  and  the  ultimate  fate  of  the  bubble  as  a  vortex  system.  Subsequent  to  this,  an 
approximate  analytical  expression  for  the  Kelvin  impulse  of  a  cavity  collapsing  in  an 
axisymmetric  geometry  near  to  a  horizontal  rigid  boundary  and  in  a  uniform 
gravitational  field  has  been  derived  by  Blake  et  al.  (1986).  The  evidence  of  the 
numerical  computations  presented  in  that  paper  indicates  that  the  direction  of  the  jet 
and  the  direction  of  the  Kelvin  impulse  at  the  end  of  the  bubble  life  are  closely 
correlated.  An  analogous  investigation  of  cavity  collapse  near  a  horizontal  free  surface 
(Blake,  Taib  &  Doherty  1987)  orovided  further  evidence  of  this  correlation. 
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In  this  paper,  the  theory  of  approximating  the  Kelvin  impulse  of  a  transient  cavity 
is  generalized  to  formally  allow  consideration  of  collapse  in  an  arbitrary  flow 
geometry.  This  theory  utilizes  the  generalization  of  the  Lagally  theorem  for  multipoles 
and  deformable  bodies,  given  by  Landweber  &  Miloh  (19X0).  It  is  shown  that  the 
hydrodynamic  force  acting  on  a  transient  cavity  is  identically  equal  to  zero,  and  this 
allows  an  expression  to  be  obtained  for  the  Bjerknes  force  exerted  upon  the  bubble  by 
the  boundaries  of  the  flow  domain  in  terms  of  the  strength  of  singularities  of  the 
velocity  potential  within  the  bubble,  defined  there  by  analytic  continuation.  Integration 
of  this  expression  over  the  approximate  lifetime  of  the  bubble  yields  an  estimate  of  the 
Kelvin  impulse. 

Particular  geometries  are  then  considered.  In  the  first  instance,  the  result  of  Blake 
et  al.  (1986)  for  the  final  Kelvin  impulse  of  a  cavity  collapsing  in  the  neighbourhood 
of  a  rigid  boundary  is  recovered.  Collapse  near  a  stationary  rigid  sphere  is  a  further 
example  considered  and  one  that  may  be  treated  exactly  in  the  approximate  theory. 
This  latter  example  introduces  the  concept  of  a  zone  of  attraction.  Motion  near  a  sphere 
in  a  uniform  gravitational  field  is  considered  and  the  final  Kelvin  impulse  estimated. 
A  surface  surrounding  the  sphere  may  be  defined  by  the  condition  that  the  line  through 
the  point  at  which  inception  of  the  cavity  occurs,  in  the  direction  of  the  final  Kelvin 
impulse,  is  tangential  to  the  sphere.  On  the  basis  of  the  proposition  that  the  directions 
of  the  final  Kelvin  impulse  and  jet  arc  closely  correlated,  this  surface  bounds  a  region 
within  which,  if  motion  occurs  inside  it.  the  jet  will  be  directed  towards  the  sphere  upon 
collapse,  and  outside  it  the  converse  is  true;  hence  the  terminology,  zone  of  attraction. 
By  utilizing  a  linearized  boundary  condition  at  a  horizontal  free  surface,  the  analysis 
is  extended  to  consider  a  sphere  submerged  to  its  centreplane.  and  the  zone  of 
attraction  is  computed.  In  the  penultimate  section  of  the  paper,  a  numerical  method  is 
described  that  allows  consideration  of  transient  cavity  collapse  in  the  neighbourhood 
of  more  general  body  shapes,  and  the  particular  class  of  shapes  considered  here  is  that 
of  bodies  of  revolution.  When  the  body  is  slender  the  results  assume  practical  relevance 
in  the  context  of  the  problem  of  underwater  explosion  bubble  collapse  near  marine 
vessels. 


2.  The  rate  of  change  of  the  Kelvin  impulse  and  the  Lagally  theorem 
The  fluid  is  considered  ideal  and  the  flow  irrotational  so  a  velocity  potential.  <P,  may 
be  introduced  that  satisfies  Laplace's  equation  in  the  flow  domain  £2.  The  boundary  of 
£2  is  P£2  and  s  S  U  E  where  5  is  the  surface  of  the  bubble  and  E  consists  of  all  other 
surfaces  that  bound  £2.  The  unit  normal  to  C£2  is  denoted  by  n  and  considered  to  be 
directed  exterior  to  £2. 

The  pressure  within  the  bubble  is  assumed  to  be  constant  throughout  the  motion  and 
equal  to  pr.  The  hydrostatic  pressure  at  the  point  of  inception  of  the  cavity  is  px.  The 
lengthscale  is  Rm,  of  the  order  of  the  maximum  bubble  radius,  and  the  timescale  is 
Rm(jj/Ap)k  with  p  the  fluid  density  and  Ap  =  p ,  -pr  the  pressure  scale.  The 
gravitational  field  is  uniform  and  given  by  g  =  -ge.,  with  e.  an  element  of  the  set  of 
Cartesian  basis  vectors.  In  non-dimensional  variables,  the  Bernoulli  equation  evaluated 
at  S  is 

~+\\V0\i  +  A-:-\  =0.  (2.1) 

<  / 

*  =  U>gRm/Ap)k 


where 

is  the  buoyancy  parameter. 


(2.2) 
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The  Kelvin 

.sc  of  a  cavity  is  (scaled  with  respect  to  /??„(/;A/>H 

/=  |  <PndS. 

J  s 

(2.3) 

and  I  changes  in 

response  to  the  action  of  forces  according  to 

dl/dt  =  F. 

(2.4) 

Integrating, 

Ht)  =  /(0)+  I  F\t')df. 

(2.5) 

and  /(O)  =  0,  so  having  detern,:»v  4  5  the  impulse  may  be  computed.  It  should  be  noted 
that  Fis  not  the  hydrodynam  *e'-  *  force  acting  on  the  bubble  due  to  the  pressure 
exerted  by  the  fluid),  which  wii,  be  Ui.ioted  by  Fh. 

The  rate  of  change  of  the  Kelvin  tpj'se  is  given  by  (Blake  &  Cerone  1982) 

dl/dt  =  F  =  F-  +  F’.  (2.6) 

with  F-  =  I  j||V0|-n-~V0jd5.  (2.7) 

and  F>  =  Ss  Ve:.  (2.8) 

F9  is  the  buoyancy  force  and  F~  the  Bjerknes  force  exerted  by  he  boundary. 
Equivalently,  the  integration  in  (2.7)  may  be  considered  as  over  the  bubble  surface  S, 
in  which  case  (Blake  &  Cerone  1982) 


Fr  =  -  j>  |^|V0|*n-^V0jdS.  (2.9) 

It  is  worthwhile  to  consider  the  relationship  between  the  above  expressions  and  the 
Lagally  theorem.  The  hydrodynamic  force  on  a  bubble  is 


and  provided  p  is  uniform  over  5.  this  is  identically  equal  to  zero!  In  the  case  of  the 
transient  cavity  considered  here,  this  is  indeed  the  case  as  the  pressure  within  the 
bubble  is  constant  throughout  the  motion.  In  view  of  this,  it  is  a  routine  matter  to 
derive  (2.6).  The  hydrodynamic  force  on  a  body,  whose  boundary  is  S,  is  given  in 
general  terms  by  (Newman  1977) 

Fh  =  -37- j>  {i|V0|*n-^V0jdS+Av-  Vez.  (2.11) 

If  this  force  is  identically  equal  to  zero  then  (2.6).  (2.8)  and  (2.9)  follow  immediately. 

The  Lagally  theorem  follows  when  the  first  two  terms  on  the  right-hand  side  of  (2. 1 1 ) 
are  evaluated  in  terms  of  the  singularities,  within  the  volume  bounded  by  5.  of  the 
potential  defined  there  by  analytic  continuation.  Landweber  &  Miloh  (1980)  have 
generalized  the  Lagally  theorem  to  treat  deformable  bodies,  and  those  for  which  the 
velocity  potential  within  the  body  has  muilipole  singularities.  In  the  course  of  their 
derivation,  Landweber  &  Miloh  (1980)  establish  the  following  expression: 

£  {i|V0|*«  -  ~  V*}  dS  =  4n  v  P(i  D°W),. 


(2.12) 
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In  deriving  this  result,  the  potential  defined  within  the  volume  bounded  by  S  is  singular 
at  points  whose  position  vectors  arc  r,  =(\\..  rx,c,).  In  the  neighbourhood  ot  the 
singularity,  the  potential  is  written  as 

0  =  0  +  *V,  (2.13) 

where  is  regular,  and  <}>  is  of  the  form 


<5«-/»,D«(I/*),  Df> 


.'‘i 


R 2  =  (  v  -  a,)2  +  (  v-y,)*  +  (r  -  ;,)*,  (2.14) 


where  the  summation  convention  has  been  used  over  q ,  the  order  of  the  multipoie,  and 
for  each  q  the  sum  is  taken  over  ait  a.  //,  y  such  that  a + (l + y  —  q.  In  (2. 1 2),  £,  denotes 
summation  over  all  singularities  within  S.  and  ( ),  denotes  evaluation  at  the  singularity 
,v;  r  s  (.r,  v,  r)  is  the  position  vector  of  an  arbitrary  point  in  space.  In  comparing  (2.12) 
with  equation  (24)  of  Landweber  &  Miloh  (1980).  note  that  the  normal  used  in  that 
paper  is  directed  oppositely  from  that  used  here. 

Using  (2.12)  in  (2.9)  gives  the  Bjerknes  force  as 


Fi  =  -4«V/>  W9i')„  (2.15) 

which  is  the  required  expression  in  terms  of  the  strength  of  the  singularities  of  <t>  within 
S,  and  derivatives  of  the  regular  part  of  the  fluid  velocity  at  these  singularities. 


3.  General  theory  for  estimating  the  Kelvin  impulse  of  a  transient  cavity 

In  what  follows  the  small  parameter  e  is  defined  to  be  0(  1  /£),  where  £  is  a  typical 
distance  from  the  bubble  centroid  to  2T,  and  the  geometry  of  L  is  considered  to  be  time 
independent.  It  is  further  assumed  that  the  buoyancy  parameter  S  is  0(e).  The  radial 
coordinate  of  the  bubble  surface  is  of  the  form 

,  i 

a(6,(j>j)  =  R(t)+'Z  ‘2  a,.„I(t)./,m(cosw)cos(",?i  +  w/.m(0).  (3.1) 

/-I  m-o 

where  is  an  associated  Legendre  function  (Gradshteyn  &  Ryzhik  1979).  This  paper 
is  concerned  with  the  lowest-order  perturbations  to  the  infinite-fluid  behaviour  caused 
by  buoyancy  and  nearby  boundaries.  The  potential  in  the  neighbourhood  of  the  bubble 
may  be  written  in  the  form  of  (2.13).  Transient  cavity  motion  in  an  infinite  fluid  has 
been  considered  by  Rayleigh  (1917),  with  the  potential  due  to  such  motion  represented 
by  a  time-dependent  source  located  at  the  centre  of  the  bubble.  This  behaviour  is 
perturbed  by  buoyancy  and  nearby  boundaries,  so  the  leading-order  contribution  to  <P 
is  a  source  term,  and  as  a  matter  of  notation  let  P„  in  (2.14)  be  written  as  m/4rc.  with 
m  the  source  strength  and  |F0|  is  0(1).  The  perturbing  influence  of  both  buoyancy  and 
nearby  boundaries  has  been  considered  by  Herring  (1950).  The  bubble  is  perturbed 
from  a  spherical  shape  throughout  its  growth  and  collapse  and  the  lowest-order 
correction  to  a  spherical  form  is  merely  a  translation  of  the  sphere  as  a  whole.  In  view 
of  this  analysis,  it  is  inferred  that  the  leading  contribution  to  |/>,|  is  0(e2),  with  the  other 
Pq  of  higher  order.  The  term  <j>’  in  (2.13)  is  non-zero  owing  to  the  presence  of  the 
boundary  L.  This  contribution  to  the  potential  may  be  considered  as  due  to  images, 
exterior  to  i2,  of  the  singularities  within  the  bubble,  and  is  0(e).  Furthermore,  the  terms 
a,  m  in  (3.1),  that  are  non-zero  owing  to  departures  from  the  spherical  shape,  are  0(e2). 
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These  mathematical  features  are  characteristic  of  the  perturbing  influence  of  both 
nearby  boundaries  and  buoyancy. 

Applying  (2.15).  the  Bjerknes  force  may  be  written  as 

F~  = -mV<y(r,.)  +  0(',ll  (3.2) 

where  rr  is  the  position  vector  of  the  bubble  centroid  and  it  is  remarked  that  the  leading 
term  in  this  expression  is  0(<  -).  The  buoyancy  force  is 

F>  =  iitrf-  /?**.  + 0(<4).  (3.3) 

since  the  terms  a,  in  (3.1)  do  not  contribute  in  evaluation  of  the  volume. 

To  proceed  further  set 

>y  =  -  ^tf(r.  rr.  s)  +  0(<r).  (3.4) 

4JI 


where  g  is  a  function  of  r.r,.  and  s  =  (.v,..v., . .v„),  which  is  a  set  of  time-independent 

distance  parameters  that  characterize  the  geometry  of  I.  Hence 

F-  =  ~  V,  ,g(  r„.  s) +  (){•'). 

4Jt 


nr 

16n 


r+tf(*'). 


(3.5) 


where  r„  is  the  position  of  the  bubble  centroid  at  t  ~  0.  Equation  (3.5)  defines  T  as  a 
function  of  r„  and  s  that  is  Note  that  f  is  defined  in  terms  of  r„.  as  rr  relates  to 
Px  (or  dipole)  contributions  to  the  potential  that  arise  owing  to  migration  of  the  bubble, 
and  since  these  are  0(k2)  the  time  variation  of  rr  contributes  to  F~  to  OUl).  As  a  further 
consequence  of  this.  rr  may  be  routinely  replaced  by  r„  in  (3.2)  and  (3.4).  and  this  is 
done  in  the  examples  considered  later  in  this  paper. 

Application  of  the  kinematic  boundary  condition  at  the  surface  of  the  bubble 
determines  the  source  strength  hi.  This  condition  is  that 


?<P  dx 
—  =  —  n. 
<  n  d  t 


(3.6) 


where  *  0/Cn  is  the  normal  fluid  velocity  and  dx/dt  is  the  velocity  of  a  Lagrangian 
surface  particle.  Using  (3.1)  to  evaluate  the  right-hand  side  of  (3.6),  and  (2.13)  to 
evaluate  the  left-hand  side,  yields 

hi  =  4nR~  R  +  0(c2),  (3.7) 

so  f- =  ir/?4/?2r+(?(f.4).  (3.8) 

In  order  to  integrate  (3.3)  and  (3.8)  with  respect  to  time  to  determine  an 
approximation  to  the  impulse,  an  expression  is  required  for  R(t).  This  is  obtained  from 
the  Bernoulli  equation,  evaluated  at  the  bubble  surface  (equation  (2.1)).  Substituting 
into  this  expression  using  (2.13)  and  (3.7)  yields 

RR  +  iRi--t-(r„.s)+  1  +0(*;2)  =  0.  (3.9) 

<  t 

Utilizing  (3.4).  and  noting  that  any  time  dependence  of  rr  contributes  to  f  1 10 
(3.9)  becomes 

RR+iRt+pRIRR  +  lW+l  =0. 


(3.10) 


so 
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with  the  error  of  OOr).  and 

!'  =  .ir(r„.4>  (3.11) 

is  #(*).  Equivalently.  (3.10)  could  be  obtained  via  an  energy  argument.  Note  that  to 
this  order  of  computation,  the  equation  giving  R(t\  is  independent  of  Multiplying 
(3.10)  by  the  integrating  factor  R-  gives 


^[hR*R-(\+,'R)\  +  R2  =  0. 


(3.12) 


At  the  maximum  radius  R  =  I ,  R  =  0.  so  integration  yields 

x  1<  1  -  /?*)*  R  H I  +,iR)  \  0  $  /  <  ir. 
-  x  1(1  -  fl'V5  R  h\+/'R)  i  $  r. 


(3.13) 


with  the  sign  of  the  square  root  of  R •  appropriately  chosen  according  to  whether  the 
bubble  is  expanding  (0  ^  t  ^  ID  or  collapsing  (\T  ^  t  $  T).  The  symmetry  of  the 
motion  about  the  maximum  radius  which  occurs  at  /  =  \T  is  noted,  with  T  the  period 
of  the  motion. 

To  proceed.  R(t)  is  determined  implicitly  and  the  Kelvin  impulse  then  evaluated. 
Now 

t  =  f  dr'  =  [H  dR’/R. 

Jo  Jo 

M. 

>R’k\+yiR)<iR’  '.T^r^T. 


^1  (\-R'3r'R’H\+l/iR')dR\  Q^t^±T. 

*•  J  0 


(3.14) 


where,  in  accord  with  the  order  of  the  approximation,  a  binomial  expansion  of 
(1  +///?')-1  has  been  utilized.  Making  the  transformation  x  -  R'3  (a  =  /?').  (3.14)  is 
written  in  terms  of  incomplete  (S^iv.r))  and  complete  (fl(ir.r))  Beta-functions 
(Abramowitz  &  Stcgun  1965) : 

(l / v  6) [B& h)  +  y>BSl. i».  0  ^  <  ir, 

/=  (2/x  6)l«M)  +  i^«i.i)J  (3.15) 

-d/\  6>l*,<«.i)+>*>(i4)J- 

At  /  =  T.  u  =  0.  so  the  period  is  given  by 

T=(2/\  6) [ff(g, i)  +  |)J.  (3.16) 

Using  (3.8),  the  Kelvin  impulse  caused  by  I  is 


=  Ttr  f  Rl  /?-  dr'  =  nr  T  R 4  RdR\ 
Jo  Jo 


(3.17) 


and  using  (3.13)  and  a  substitution  as  above  to  carry  out  the  integration  yields 

A v  6rt/9) r[fi,(|, i)-ipU,(i. |)],  ^  IT, 

/-(/)={( 2\  6n/9) j)  —  \/iB(t |)J  (3.18) 

- ( x  6nl9) r[B,(l  1) - &BJLI §)],  \T^  t 
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wtth  an  error  of  OU.'l.  A  similar  evaluation  yields  the  impulse  due  to  buoyancy  as 
C\  6ic/S»A'- +  0  $  /  ^  \T. 

/"(/>  =  (4X  6n/9).la[«V.i)  +  ^«T.iHr-  (3.19) 

-(2V  kT^i^T. 

with  an  error  of  0(*‘). 

In  order  to  estimate  the  period  and  Kelvin  impulse  of  a  transient  cavity,  the 
coefficients  p  and  T  must  be  evaluated.  These  are  purely  functions  of  the  geometry  of 
I  and  thus  characterize  its  influence  upon  the  motion  of  the  cavity.  Several  examples 
in  the  next  sections  illustrate  the  approach. 


4.  Examples  treated  exactly  in  the  approximate  theory 

4. 1 .  ( 'olhipse  near  a  horizontal  rigid  plane 

Consider  transient  cavity  collapse  near  a  horizontal  rigid  plane  at  r  =  0.  The  regular 
potential  in  the  neighbourhood  of  the  bubble,  to  OU.).  is  due  to  an  image  source  of 
strength  m  reflected  about  the  plane.  Considering  inception  of  the  cavity  to  occur  at 
r„  s  (0, 0.  £). 


m 


<!>’(')  =  -4^<r'r»>‘ 


(4.1) 


with 


X(rr,,)  = 


lr -(0.0. -01‘ 


From  (4.2)  F  and  p  are  evaluated  as 


(4.2) 


r=--e. 


(4.3) 


where  e<  is  a  unit  vector  directed  away  from  the  boundary.  Blake  el  at.  ( 1986)  did  not 
consider  the  correction  to  the  Rayleigh  equation  for  R  due  to  the  nearby  boundary,  so 
to  recover  the  result  of  that  paper  p  is  set  equal  to  zero  in  (3.18)  and  (3.19). 

Since  the  Kelvin  impulse  is  the  impulse  that  would  have  to  be  applied  over  the 
surface  of  the  bubble  in  order  to  generate  the  observed  flow  Held  from  rest,  there  is 
some  physical  basis  for  supposing  that  the  direction  of  the  jet  and  the  direction  of  the 
Kelvin  impulse  at  the  time  of  collapse  are  closely  correlated.  The  analytical  estimates 
of  the  impulse  at  the  end  of  the  bubble  collapse  obtained  by  Blake  et  al.  (1986.  1987) 
were  found  to  provide  reasonable  estimates  of  the  Kelvin  impulse  of  the  jet-pierced 
bubble,  despite  the  analysis  assuming  small  departures  from  a  spherical  shape.  This 
follows  from  the  speed  with  which  jet  formation  occurs.  Typically,  the  time  over  which 
this  non-spherical  collapse  occurs  is  of  the  order  of  2  3%  of  the  bubble  life.  Over  some 
interval.  At.  the  change  in  the  Kelvin  impulse  is 


A  / 


rt+At 

=  l  1 


FU‘)dt\ 


(4.4) 


so  that 


jAij  ^  At  sup  1/^01, 


(4.5) 


where  sup  denotes  the  supremum  value.  F  is  bounded  throughout  the  motion  so  that, 
over  a  sufficiently  small  time  interval,  the  magnitude  of  the  impulse  change  is  small. 
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The  time  at  which  the  bubble  enters  the  rapid  collapse  phase  and  becomes  highly  non- 
spherical  is  given  with  good  accuracy  by  (.V16).  and  up  until  this  time  the  Kelvin 
impulse  can  be  determined  using  the  results  of  §  3  Because  of  the  generally  short  time. 
At.  over  which  the  collapse  occurs,  this  provides  a  reliable  estimate  of  the  impulse  of 
the  non-spherica!  collapsed  bubble,  the  change  in  the  impulse  brought  about  by  the 
deformation  from  a  spherical  shape  being  small,  in  accord  with  (4.5). 


4.2.  Collapse  near  a  stationary  rigid  sphere 

Without  loss  of  generality,  suppose  that  the  centre  of  a  fixed  rigid  sphere  of  radius  a 
is  located  at  the  origin  and  that  inception  occurs  at  r  =  £.  In  this  case  is  given  by  the 
Weiss  sphere  theorem  (see  Milne-Thomson  I960) 


am  I  m  £  dr 

4n£  \r  -  (0. 0.  a1/ £)|  +  4iui  J„  |r-(().().7)f 


(4.6) 


The  first  term  corresponds  to  a  source  of  strength  <//»/£  located  at  r  =  </-/£,  and  the 
second  term  to  a  uniform  linear  distribution  of  sinks  of  density  m/a.  from  r  =  0  to 
r  =  «*/£.  This  is  the  image  set  associated  with  a  source  located  at  r  =  £  and  is  such  that 
the  boundary  condition  V&-n  =  0  is  satisfied  at  the  surface  of  the  sphere.  In  this 
example  the  geometry  of  I  is  characterized  by  s  =  a.  From  (4.6) 


a  I  I  |V£  dr 

£|r  — (0.0.  <r/£)|  </J„  |r  — (O.O.r)l’ 


(4.7) 


and 


r=- 


4«a 


(4.8) 


where  the  onl>  non-zero  component  of  the  Bjerknes  force  exerted  by  the  sphere  is 
directed  towards  the  origin,  as  denoted  by  use  of  the  radial  vector  er. 

Since  the  direction  of  the  Kelvin  impulse  at  the  conclusion  of  the  collapse  is  expected 
to  be  closely  correlated  with  the  direction  of  the  jet.  the  theory  developed  here  may  be 
used  to  infer  the  behaviour  of  buoyant  cavities  collapsing  in  the  neighbourhood  of  a 
sphere.  The  angle,  i//.  that  the  direction  of  the  final  impulse  makes  with  the  line  joining 
the  point  of  inception  to  the  centre  of  the  sphere  is 

$  =  arctan  (/„(  Tr)/Ir{  Tr)).  (4.9) 

where  /„  and  Ir  are  the  components  of  the  Kelvin  impulse  in  spherical  coordinates.  The 
geometry  is  shown  in  figure  1 .  The  critical  angle.  yr.  is  defined  by 

i/>r  =  arcsin  («/£).  (4.10) 

and  t  =  Vc  (4.11) 

defines  a  condition  such  that  the  direction  of  the  final  Kelvin  impulse  is  tangential  to 
the  sphere.  Since  /( Tr)  is  a  function  of  £  and  0.  (4.11)  defines  a  surface  surrounding  the 
sphere  and  this  surface  shall  be  hereafter  denoted  by  S„.  Within  this  surface  the 
Bjerknes  attraction  of  the  sphere  is  sufficiently  large  that  the  final  Kelvin  impulse  is 
directed  towards  it.  It  is  supposed  that  the  jet  formed  will  be  similarly  directed.  Thus 

the  terminology  :one  of  attraction  is  introduced  to  describe  this  region  of  space. 

Conversely,  outside  the  surface  defined  by  (4. 1 1 ).  the  final  Kelvin  impulse  is  directed 
away  from  the  sphere. 

The  example  a  =  2.0,  8  =  0. 1  is  shown  in  figure  2.  Equation  (4. 1 1 )  has  been  solved 
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Point  ol'  inception 

FiGURH  1 .  The  geometry  for  buoyant  cavity  motion  near  a  stationary,  rigid  sphere. 


Figure  2.  The  zone  of  attraction  around  a  sphere  of  radius  a  =  2.0. 

The  buoyancy  parameter  is  =  0. 1 . 

by  the  bisection  method,  with  substitution  of  (4.8)  into  (3.18)  and  (3.19)  yielding  the 
appropriate  components  of  the  impulse.  Sa  encloses  the  sphere  and  as  r  -» —  oo  the 
section  of  this  surface  in  some  plane  z  =  constant  tends  to  xl  +  yi  =  a1.  This  is  expected 
on  physical  grounds.  As  z  -*  —  oo  the  presence  of  the  sphere  is  little  felt  {T ,/r  -*■  0)  so  the 
dominant  buoyancy  force  yields  a  final  Kelvin  impulse  directed  upwards  and  towards 
the  sphere,  if  the  point  of  inception  is  under  it. 

In  regard  to  jet  formation  as  a  damage-causing  mechanism,  the  zone  of  attraction 
is  relevant  as  it  defines  a  region  in  which  jet  impact  may  occur.  However,  the  part  of 
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this  zone  extending  to  :  =  -  /  is  of  little  relevance  in  this  context  as  a  jet  forming 
towards  the  sphere,  but  a  long  way  from  it.  will  cause  no  damage. 

4.3.  Collapse  near  a  stationary  sphere,  submerged  to  its  centreplane  at  a  free  surface 

Consider  now  a  sphere  submerged  to  its  centreplane  at  a  horizontal  free  surface  defined 
by  r  =  0.  To  the  accuracy  of  computation  required,  the  lineartzed  boundary  condition 
that  the  potential  vanishes  on  the  free  surface  is  applied.  For  the  purpose  of  calculation 
of  r.  axisymmetry  may  be  exploited  and  *he  point  of  inception  considered  as 

r„  =  (w.0.f).  (4.12) 

with  the  centre  of  the  sphere  the  origin.  The  distance  between  the  point  of  inception 
and  the  sphere  centre  is 

£  =  (,.,- +  £*)->.  (4.13) 

The  simplest  way  to  satisfy  the  boundary  condition  on  :  =  0  is  to  introduce  an  image 
sink  of  strength  m  at 

r,  =  (w.0.  -£).  (4.14) 

which  is  the  reflection  of  r„  about  the  plane  r  =  0,  and  determine  0  such  that 
V<t>n  =  0  on  the  whole  of  the  sphere  surface.  From  symmetry  considerations,  it  is 
evident  that  the  resultant  potential  will  satisfy  0  =  0  on  :  =  0. 

The  sphere  theorem  may  be  applied  to  both  the  source  representing  the  bubble  and 
its  image  sink  to  give 


$\r)  = 


4rt|r  — (w.0,  -£)| 


r  and  //  arc 


am  I  m  f"' £  dr 

“ 4*£ |r-((iv,>/f-.().<rC/£-)|  +  4™  J „  |r -(.■„,/£.  0.  r£/£) 

am  l  m  |V£  dr 

+  4^£  I r  0-  -  «W)|  _4™J„  |r  — (rw/£.  0,  -r£/£)f 


r=  l\  er  +  r.  e.. 

4 a'o>  2  j 4a-,.rf-  -(</-  +  £-)(£'  +  <<•  +  2<r£s) 


r  =  ~?(?-aX  +  m,\  (^  +  u‘  +  2u*(r-or))*  . 


(4.15) 


(4.16) 

(4.17) 


^  4 aX  ,  2  \4aXV-(<r -?)&  +  «' ~2aior)  ,].sgn(£)  ,4  104 

l  crtffW))* - 'r-^-  (4,8) 

1  ^  «  i,  i  e  i 

11  2|£|  +  f-,r  Jn|£W/ 


(£4  +  rt4  +  2a-(£2-w-))!  +  «ln 


1 .  l(£'  +  «'  +  2a-(£---^-))4r/-  +  r--^) 


(4.19) 


where  er  and  e,  are  cylindrical  basis  vectors. 

The  zone  of  attraction  is  defined  by  the  condition  that  the  line  passing  through  the 
point  of  inception,  in  the  direction  of  the  final  Kelvin  impulse,  is  tangential  to  the 
surface  of  the  sphere.  The  equation  of  this  line  is 

r'  =  o,+  \Ir{Tr ).  =  t  +  AI:(Tr). 


(4.20) 


Figure  3.  The  zone  of  attraction  around  a  sphere  submerged  to  its  centrcplane  at  a  free  surface.  The 
cross-section  through  the  axis  of  symmetry  is  shown,  rf  =  0.1.  The  radius  of  the  sphere  is  («)</=  1 .0, 
(b)  a  =  1.2.  (<  )  a  =  1.4.  (</)  a  =  2.0. 

where  /\  is  a  parameter  and  the  components  of  the  linal  Kelvin  impulse  arc  functions 
of  o>  and  C,  which  detine  the  point  of  inception.  The  equation  detining  the  surface  of 
the  sphere  is 

(4.21) 

so  that  the  points  at  which  the  line  and  sphere  intersect  are  given  by  the  solutions  of 

(/;(rf)  +  /;(rr))A-  +  2Mr(^  +  a.(rr)),\+(f -<r)  =  0.  (4.22) 

For  the  line  to  be  tangential  to  the  sphere  this  equation  has  only  one  solution,  so  the 
equation  defining  the  surface.  S„,  that  bounds  the  zone  of  attraction,  is 

Aok  Q  =  («/,(  T; >  +  £/,<  Tr)f  -  ( r;(  Tr)  +  /r(  Tr))(£-  - <r)  -  0.  (4.23) 

This  function  is  negative  outside  the  zone  of  attraction  and  positive  inside  it.  Utilizing 
this  fact,  (4.23)  has  been  here  solved  using  the  bisection  method.  Note  also  that  for  the 
solution  to  be  physically  relevant  A  must  be  positive,  indicating  that  the  impulse  is 
directed  towards  the  sphere. 

Examples  are  shown  for  S  =  0.1  and  various  values  of  a  in  figure  3.  The  section 
through  Sa  in  the  axis  of  symmetry  is  shown  for  a  —  1.0  (figure  3 a),  a  =  1.2  (figure  3  b), 
a  =  1.4  (figure  3c)  and  a  =  2.0  (figure  3 cl).  Consider  first  a  =  1.0.  The  striking  feature 
is  that  the  zone  of  attraction  consists  of  two  distinct  regions.  This  may  be  explained  by 
considering  the  forces  acting  on  the  bubble.  In  close  proximity  to  the  sphere  its 
Bjerknes  attraction  is  dominant  causing  the  jet  to  be  directed  towards  it,  noting  that 
the  effect  of  the  free  surface  is  to  cause  jet  formation  downwards.  As  the  point  of 
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inception  moves  downwards,  the  Bjerknes  force  exerted  by  the  sphere  falls,  and  the 
effect  of  the  free  surface  becomes  dominant  and  the  final  Kelvin  impulse  is  directed 
away  from  the  sphere.  At  greater  depths  the  free  surface  and  Bjerknes  force  exerted  by 
the  sphere  are  little  felt,  and  buoyancy  causes  jet  formation  upwards  and  towards  the 
sphere.  As  r  -  -  x  the  section  of  Sn  in  the  plane  r  =  constant  tends  to  that  of  the 
sphere  in  the  plane  z  =  0. 

A  further  significant  feature  is  that  the  bottom  section  of  the  zone  exhibits  a  branch 
that  tends  outwards  and  is  asymptotic  to  the  plane  defined  by 


2<J2 


=  0. 


(4.24) 


The  significance  of  this  plane  is  that  in  the  absence  of  the  sphere,  it  defines  that  depth 
of  inception  at  which  the  final  vertical  component  of  the  Kelvin  impulse  is  zero  in  the 
approximate  theory.  In  this  case,  for  <o  sufficiently  large,  just  below  the  surface  defined 
by  (4.24)  there  is  a  small  upwards  component  of  the  impulse  that  combines  with  the 
small  component  parallel  to  the  free  surface  and  due  to  the  presence  of  the  sphere,  in 
order  that  the  final  impulse  is  directed  towards  the  sphere.  Note  also  that  the  part  of 
Sa  near  to  the  sphere  is  only  calculated  to  some  close  point,  as  the  theory  is  not 
applicable  in  the  very  close  neighbourhood  of  its  surface. 

The  changing  structure  of  the  zone  of  attraction  as  a  increases  is  interesting  and  the 
sequence  shown  in  figure  3  illustrates  this.  For  a  -  1.2  the  larger  sphere  causes  a  greater 
Bjerknes  force  so  that  part  of  the  zone  of  attraction  immediately  below  the  sphere 
extends  to  a  greater  depth.  Figure  3(c)  (a  =  1.4)  illustrates  the  merging  of  the  upper 
and  lower  portions  of  the  zone  of  attraction.  In  this  case  the  sphere  is  sufficiently  large 
that,  near  to  the  r-axis  at  least,  there  is  no  region  where  the  repulsion  of  the  free  surface 
is  dominant.  The  example  of  a  —  2.0  (figure  3c/)  shows  that  for  a  sphere  reasonably 
larger  than  the  bubble,  the  free  surface  is  significant  only  at  the  edges  of  the  sphere. 


5.  The  numerical  method 

To  exploit  the  full  power  of  the  theory  developed  here  for  estimating  the  Kelvin 
impulse  of  a  transient  cavity,  it  is  necessary  to  investigate  numerical  methodologies  for 
dealing  with  complicated  flow  field  geometries.  In  order  to  evaluate  the  approximation 
to  the  Kelvin  impulse,  suppose  that  the  regular  part  of  the  potential  in  the 
neighbourhood  of  the  bubble  centroid  is  considered  as  due  to  a  distribution  of  sources 
of  surface  density  mcr,  over  Z.  Then 


=  -^g(r,ru,j) 

(5.1) 

with 

(5.2) 

Hence 

r.  4  f 

Jz  \r0-<l\ 

(5.3) 

t- 

and 

f 

J  z  K  91 

(5.4) 

Note  that  a  depends  upon  the  position  on  Z  (the  vector  9),  the  vector  s  that 
characterizes  the  geometry  of  Z,  and  the  point  of  inception  r„. 
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Writing 


0  =  - 


m _ l_ 

4rt  I  r  -  r,. 


and  restricting  attention  to  cases  where  £  is  rigid,  the  boundary  condition 

V<tHp)n(p)  =  0.  (5.6) 

satisfied  for  all  pel.  determines  fr{q.r„,s).  Using  (5.5)  in  (5.6),  the  following  integral 
equation  is  obtained  for  rr : 

C  (5.7) 

IP-'-.l  Jl  Ip-fl 

A  particular  difficulty  with  this  equation  is  the  singularity  in  the  integrand  at  q  —  p\ 
however,  this  can  be  removed  (Landweber  &  Shahshahan  1992)  by  noting  that 


Then  (5.7)  may  be  written  as 

,n  _  r>  _  1  ip-ra)-n(p) 
n  \p-rf 


If  M<l^0.s)[p-q)-n{p)-fT{p.r0 
4nJ.  \p-q\3  " 


$)(?-/>)•«<?)} 


d5(9).  (5.9) 


Equation  (5.9)  is  a  Fredholm  integral  equation  of  the  second  kind  and  provided  the 
process  converges,  may  be  solved  by  iteration  as  follows: 


<r{m(p.ri„s)  =  — 


1  (/>-/•„)•«(/») 


4rt  \p-r„\ 


(T{n)(p.r0.s)  =  <rm(p,  r„.  s) 

1  f  r„,  s){p-q)- njp)  -rrln-,){p.  r„,s)(q  -p)  •  n{q)) 

4njr  \p-q\3 


(5.10) 


dS(q).  (5.11) 


where  the  superscript  (»)  denotes  the  approximation  obtained  by  performing  n 
iterations  of  (5.10)  and  (5.11).  The  particular  numerical  details  in  carrying  out  this 
iteration  may  vary  depending  upon  the  particular  geometry. 

To  illustrate  one  approach,  the  specific  example  where  E  is  a  body  of  revolution 
about  the  .v-axis  will  be  considered  in  detail.  The  flow  domain  is  infinite  and  the 
position  of  the  body  fixed.  The  surface  of  the  body  is  given  in  parametric  form  as  a 
function  of  £  and  6  as 

p  =  (A*(£),T(£)sin0,T(£)cos0),  (5.12) 

where  6  is  measured  from  the  vertical.  The  variable  £  is  the  arclength  along  the  curve 
that  generates  the  body.  The  radius  of  the  body,  r,  may  be  equivalently  considered  as 
a  function  of  X .  in  which  case  it  is  evident  that 


£  =  J*V(1  +r*v(A"))^dAr', 


(5.13) 
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Figure  4.  The  geometry  of  a  body  of  revolution. 


the  subscript  .V  denoting  differentiation.  The  geometry  is  shown  in  figure  4  and  as  a 
consequence  of  the  parameterization  with  respect  to  £. 

X;  +  ti  =  I.  (5.14) 

Axisymmetry  may  be  exploited  here  in  the  evaluation  of  r  and  //.  so  inception  is 
considered  to  occur  at 

r„  =  (£,0.w).  (5.15) 

The  unit  normal  to  2T,  exterior  to  the  flow  domain  is  given  by 

n  »  (r4,  —  Xc  sin  0,  —  X<  cos  V).  (5.16) 

and  the  surface  area  element  is 

dS  =  rd0d£.  (5.17) 

Using  (5.12)  and  (5.15H5.17)  in  (5.9)  yields 

*s.0)  =  1 -ri<oco.<»)-fri(0(TO-o-rto^o+ j_  p-*d{. 

( W£)  —  0*  + «-  -  2wt(£)  cos  (0)  +  r-(f))'- 
x  [*•(£',  O')  {( Af(£)  -  A '(£'))  Tjd)  -  *>£)  r(£)  +  ,Y£(£)  r(£')  cos  (0 - 0')} 

-  o(£,  0)  {(*(£')  -  X(0)  -  A>(£')  r(f)  +  *,<£')  r(£)  cos  (0  -  0')}]/ 

[WO  -  *(£'))*  +  r*(£)  +  T-(0  -  2  r(£)  t(£')cos  (0-0)]*.  (5.18) 

Symmetry  about  the  plane  y  =  0  is  exploited  in  order  to  reduce  the  interval  of  0' 
integration  to  [0,n].  The  interval  [£mi„,£mnJ  is  partitioned  into  /if  segments  by  the 

sequence  {£,,£* . £n{+,}.  Similarly,  [O.n]  is  partitioned  into  n„  segments  by  the 

sequence  {0,.0, . 0n?+1}.  Notice  that  £,  =  £min,  -  £mox,  0,  =  0  and  0^,  =  rt. 

The  following  quantities  are  also  defined: 


^£t+i  —  £<+i  £<»  * —  K2. . /if.  | 

^+I=0;+|-0;,  >=1.2 . rt„.) 


(5.19) 


This  division  defines  a  mesh  over  £moxJ  x  [0,  it]  which  consists  of  tun,,  rectangular 
elements.  Element  ij  has  vertices  (£„0;),  (£,+I,0;).  (£.>,,0;.,).  (£*,0;.,).' 

The  iteration  of  (5. 18)  following  (5. 10)  and  (5. 1 1 )  is  carried  out  at  a  set  of  ncn„  nodes 
defined  by 

(£m^h)  =  (£<+^£<+„^+^+1),  /=l.2 . /»,.  >=1.2 . n„.  (5.20) 

In  order  to  perform  the  integration  over  element  /.>  in  (5.18),  the  source  density  <t  is 
assumed  constant  and  equal  to  <ru ,  its  value  at  (£(+j,0;+p.  A  product  Gauss  rule  is 


Tlw  Kelvin  impulse  of  a  transient  cavity  89 

employed  with  an  even  number  of  abscissae,  in  order  that  no  Gaussian  point  coincides 
with  a  node  point,  in  order  to  evaluate  r.  .V.  t£  and  .V£  at  the  Gaussian  points,  a  cubic 
spline  is  used  to  approximate  t(£)  and  A'(£)  and  constrained  to  pass  through  the  values 
of  these  functions  at  £,  \i  =  1.2 . //£+l|.  Evaluating  (5.18)  utilizing  these  rep¬ 

resentations  yields  the  equation 


••  i  a  „ 

_1(>I  .  V  V  n  rri"~u 
"a  "a  T  —I  '■'an  fTki 


(5.21) 


Gljkl  ~  A.  -j  "t.i, 

Kj^l) - A'(£,  J) t£(£, - ,V£(£, .;) r(£, f j) 4- ■(',(£,.;) r(£t.  J cos +  ( -  1  f  <>K,)\ 
((*<&♦})  -  •*'(£*.  ,,))*  +  r2(& , . )  +  t'(£,. -  2r(&  ,*)  r(^.. ,{)  cos  (0, .  j  +  ( -  I V <„))’ 

I  i  ",  m-,  /<„„ 

-JL.fi  rf  V  V  V  V  V  I.-  Tic  \ 

JZnlk  °ll  —I  —  —I  -d  nk.i.nU.Tl*k.i.l 
j-l  k-l  /-I  ij-l  i.-l 

[(AT(^,J-A'(^+,))r£(^J-.V£(^,Jr(^,J  +  A^^,Jr(f,,)cos(«>t.+(-l)^,.J] 


[(*<£,*}) -  AT(|,  „))a  +  r-(^+i)  +  r-(|,. „) - 2r(£t>1) r(£*, „)cos  (tf, 4-  ( -  1  )* 0,. ,„)]’ 

(5.22) 

In  this  expression.  ngi  and  ng„  respectively  denote  the  number  of  Gaussian  points  used 
to  carry  out  the  integration  over  £'  and  0\  wk.  it  and  w,  im  denote  the  weights  and  £k  i(  and 
f),  iit  the  abscissae  for  integration  over  the  element  k.  I.  Slk  is  the  Kronecker  delta,  as  is 

To  test  and  illustrate  the  numerical  method,  consider  the  oblate  spheroid  given  by 
revolution  of  the  curve 

A'*  +  r7//-  =  (5.23) 

about  the  .v-axis.  with  //  ^  1 .  The  case  //  —  I  is  a  sphere  and  comparison  with  the  exact 
solution  presented  in  the  previous  section  allows  the  numerical  scheme  to  be  validated. 
In  order  to  partition  the  interval  [£m(fl,£ma,J.  points  are  chosen  as 

X,  =  /  sin (*'/,)•  rf  =  «/ cos(*7).  |  ^5  24) 

i/r,  =  -hn  +  (i- \)n/nt,  /=  1,2 . /ij+l.j 

A  cubic  spline  is  fitted  to  these  points  using  the  technique  of  Kucera  (1993)  that  yields 
the  parameterization  with  respect  to  the  exact  arclength  £.  A  particular  benefit  of 
utilizing  (5.24)  to  partition  [£mi„,ZmaA  is  that,  in  the  case  where  the  body  is  slender, 
smaller  elements  are  chosen  near  the  ends  where  the  curvature  is  high.  In  this  example 
M  is  a  constant. 

The  performance  of  the  numerical  method  in  the  case  of  a  sphere  of  radius  1  is 
illustrated  by  figure  5 (a,  b).  Figure  5(a)  shows  the  relative  error  in  Tn  the  non-zero 
component  of  T  directed  towards  the  centre  of  the  sphere,  and  figure  5(b)  shows  the 
relative  error  in  //.  The  variation  of  the  error  with  £.  the  distance  of  the  inception  point 
from  the  centre  of  the  sphere,  is  shown  for  various  «£  x  n„  combinations,  in  particular 
14x14,  19  x  19.  24  x  24.  29  x  29  and  34  x  34.  In  all  computations  ngl  =  ng„  =  2  and 
(5.21)  was  iterated  until  the  solution  had  converged  in  the  sixth  decimal  place.  This 
required  approximately  7  to  1 5  iterations,  the  larger  number  when  £  *  1 .  These  graphs 
display  the  convergence  of  the  method  as  the  mesh  is  refined.  Furthermore,  they 
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Figure  5.  The  relative  error  in  calculation  of  (a)  f,  and  (b)  /<  for  the  example  of  a  sphere  of  radius 
I.  The  variation  of  this  error  with  £  is  shown  for  the  mesh  sizes:  (i)  14  x  14.  (ii)  19  x  19.  (iii)  24  x  24. 
(iv)  29  x  29.  (v)  34  x  34. 

illustrate  the  divergence  as  £-*•!.  This  corresponds  to  the  point  of  inception 
approaching  the  sphere  surface  and  in  view  of  this  interpretation  it  is  evident  that  the 
numerical  calculation  should  diverge  under  these  circumstances.  If  calculations  were 
required  to  very  high  accuracy  in  the  event  of  £  as  1 .  the  method  could  be  improved  by 
noting  that  in  such  a  case  the  flow  field  in  the  neighbourhood  of  the  bubble  source 
appears  as  for  a  source  near  a  plane  boundary,  an  example  for  which  the  solution  is 
known.  Alternatively,  a  mesh  refinement  process  could  be  utilized.  However,  the 
validity  of  the  theory  is  questionable  in  such  a  regime  so  it  does  not  seem  appropriate 
to  invest  a  great  deal  of  effort  in  highly  accurate  calculations  there. 

The  numerical  method  for  calculation  of  f  and  n  may  be  used  to  compute  the  zone 
of  attraction  around  various  bodies,  just  as  was  done  for  the  example  of  a  sphere.  The 
zone  of  attraction  around  the  sphere  shown  in  figure  2  was  calculated  using  the 
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Figure  6.  The  zone  of  attraction  around  the  body  of  revolution  defined  by  (5.23).  /  =*  10.  //  =  0.2. 

<5  »  0.1.  The  cross-sections  through  the  planes  {a)  x  =  0  and  (/>)  y  =  0  are  shown. 

numerical  method  employing  nt  =  n„=  14  and  10  iterations  of  (5.21).  The  result  was 
visually  indistinguishable  from  that  computed  using  the  analytical  expressions.  Since 
the  smallest  value  Ha  takes  is  about  2.  inspection  of  figure  5 (a.  b)  reveals  that  both  Tr 
and  p  are  computed  with  a  relative  error  of  less  than  I0‘2.  and  this  has  proved  to  be 
entirely  sufficient. 

As  an  example  relevant  to  the  problem  of  underwater  explosion  bubble  collapse  near 
a  marine  vessel,  consider  the  zone  of  attraction  about  the  body  defined  by  (5.23)  in  an 
infinite  fluid,  with  a  uniform  gravitational  field  acting.  The  geometry  of  Sa,  in  the 
planes  .v  =  0  and  y  =  0,  is  shown  in  figure  6  for  the  case  f  —  10. »/  =*  0.2  and  8  —  0.1. 
In  this  example  n,  =  19,  and  n„  =  14  with  ten  iterations  of  (5.21)  carried  out.  As  for 
collapse  near  a  sphere,  as  r  —  oo  the  section  of  S„  in  the  plane  r  =  constant  tends  to 
the  section  of  the  body  in  the  plane  z  =  0. 

In  concluding  this  section,  it  is  noted  that  the  CPU  time  required  to  compute  an 
approximation  to  the  Kelvin  impulse  using  the  technique  described  here  is  less  than  the 
time  taken  to  perform  one  time  iteration  of  a  boundary  integral  calculation  of  the 
whole  flow  field  (see  for  example  Harris  1992).  This  provides  an  illustration  of  the 
power  of  this  theory.  An  estimate  of  the  important  parameter  of  the  direction  of  jet 
formation  upon  transient  cavity  collapse  may  be  made  with  considerable  economy  in 
some  rather  complex  but  practically  relevant  geometries. 


6.  Conclusions 

In  this  paper,  an  expression  has  been  given  that  relates  the  components  of  the 
Bjerknes  force  acting  on  a  bubble  to  the  singularities  of  the  fluid  velocity  potential 
within  the  bubble,  defined  there  by  analytic  continuation.  Application  of  this 
expression  to  transient  cavities  allows  a  general  theory  to  be  developed  for  estimating 
the  Kelvin  impulse  during  collapse  in  a  flow  field  of  general  geometry,  assuming  that 
perturbations  from  spherical  shape  are  small.  Since  the  large  deformation  from  a 
spherical  shape  brought  about  by  jet  formation  occurs  quickly,  the  estimates  obtained 
assuming  small  departures  from  a  spherical  shape  are  expected  to  well  approximate  the 
impulse  of  the  cavity  after  jet  formation  has  occurred.  Using  this  theory,  the  result  of 
Blake  et  al.  (1986)  for  the  impulse  of  a  cavity  collapsing  near  a  plane  rigid  boundary 
has  been  recovered,  and  the  effort  required  in  carrying  out  the  analysts  is  significantly 
less  than  for  the  original  approach. 


92 


J  P.  Best  and  J.  R.  Blake 

The  power  of  ihe  general  theory  has  been  illustrated  by  consideration  of  the  further 
examples  of  transient  cavity  collapse  near  a  sphere,  both  in  an  infinite  fluid  and 
partially  submerged  at  a  horizontal  free  surface.  Furthermore,  a  numerical  im¬ 
plementation  of  the  method  has  been  illustrated  that  allows  complex  flow  field 
geometries  to  be  considered.  Since  the  direction  of  jet  formation  upon  the  collapse  of 
a  transient  cavity  appears  to  be  correlated  with  the  direction  of  the  Kelvin  impulse,  the 
analytical  and  numerical  estimates  obtained  for  the  impulse  may  be  used  to  predict  the 
direction  of  jet  formation.  In  this  context  the  concept  of  a  zone  of  attraction 
surrounding  an  isolated  body  has  been  introduced,  if  inception  of  a  cavity  occurs  in 
this  region  of  the  flow  domain,  the  Anal  Kelvin  impulse  (and  presumably  the  jet)  will 
be  directed  towards  the  body.  The  bounding  surface  of  this  region  may  be  calculated 
using  the  estimates  of  the  final  Kelvin  impulse  computed  in  this  paper. 

Significant  assumptions  have  been  made  in  developing  the  theory  of  this  paper,  as 
stated  above.  However,  the  evidence  obtained  to  date,  particularly  from  numerical 
simulations,  indicates  that  the  results  derived  here  provide  a  reliable  indicator  of  jetting 
behaviour  upon  the  collapse  of  transient  cavities.  The  numerical  simulations  of  Blake 
et  al.  (1986.  1987)  indicate  well  the  qualitative  correlation  between  the  directions  of  the 
Kelvin  impulse  at  collapse  and  the  jet.  These  results  further  indicate  that  analytical 
approximations  derived  assuming  that  t  is  small  still  provide  reasonable  predictions  of 
jetting  behaviour  when  inception  of  a  transient  cavity  occurs  at  approximately  one 
maximum  bubble  radius  from  a  rigid  boundary  or.  in  the  notation  of  this  paper,  when 
ft  is  0(  I ).  As  a  further  comment,  cavity  collapse  near  a  compliant  boundary  has  been 
considered  briefly  by  Blake  (1988).  and  the  analytical  estimate  of  the  final  Kelvin 
impulse  provides  predictions  of  jetting  behaviour  that  are  in  accord  with  the  scarce 
experimental  data  available. 

In  order  to  determine  fully  the  realms  of  applicability  of  the  theory  developed  here, 
further  experimental  and  numerical  work  is  required.  Both  are  currently  proceeding, 
especially  addressing  cavity  collapse  in  the  neighbourhood  of  the  non-trivial  geometries 
considered  here. 

We  would  like  to  thank  an  anonymous  referee  whose  constructive  comments  have 
contributed  to  a  substantially  improved  Anal  manuscript. 
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